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Abstract: High-contrast imaging of extrasolar planet candidates around a 
main-sequence star has recently been realized from the ground using current 
adaptive optics (AO) systems. Advancing such observations will be a task 
for the Gemini Planet Imager, an upcoming "extreme" AO instrument. 
High-order "tweeter" and low-order "woofer" deformable mirrors (DMs) 
will supply a >90%-Strehl correction, a specialized coronagraph will sup- 
press the stellar flux, and any planets can then be imaged in the "dark hole" 
region. Residual wavefront error scatters light into the DM-controlled dark 
hole, making planets difficult to image above the noise. It is crucial in this 
regard that the high-density tweeter, a micro-electrical mechanical systems 
(MEMS) DM, have sufficient stroke to deform to the shapes required 
by atmospheric turbulence. Laboratory experiments were conducted to 
determine the rate and circumstance of saturation, i.e. stroke insufficiency. 
A 1024-actuator 1.5-/im-stroke MEMS device was empirically tested with 
software Kolmogorov-turbulence screens of ro =10-15 cm. The MEMS 
when solitary suffered saturation ~4% of the time. Simulating a woofer DM 
with ~5-10 actuators across a 5-m primary mitigated MEMS saturation 
occurrence to a fraction of a percent. While no adjacent actuators were 
saturated at opposing positions, mid-to-high-spatial-frequency stroke did 
saturate more frequently than expected, implying that correlations through 
the influence functions are important. Analytical models underpredict the 
stroke requirements, so empirical studies are important. 
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1. Introduction 

Direct imaging of extrasolar planet candidates around main-sequence stars has recently been 
realizedim |2] [51 employing current adaptive optics (AO) systems and the angular differential 
imaging technique JU. Such high-contrast observations with ground-based telescopes require 
careful control of the incident wavefront, eliminating noise so that faint sources can be detected. 
Designed specifically for high-contrast imaging, the Gemini Planet Imager (GPI)0|6] will use 
an "extreme" adaptive optics (ExAO) system to flatten the incident wavefront. 

Key to this ExAO system is the high-order deformable mirror (DM), a polysilicon 
micro-electrical mechanical systems (MEMS) DM developed by Boston Micromachines 
Corporation f?!. Consisting of thousands of individually addressable actuators in a coin-sized 
space, MEMS DMs are a low-cost approach to high-order wavefront correction; however, 
MEMS devices exhibit less mechanical stroke than conventional piezo-actuated DMs. 

Stroke saturates when the desired DM shape at a given actuator is beyond that actuator's 
range of motion. Saturation degrades the wavefront correction and scatters light into the region 
that should remain dark for detection of faint sources. This paper evaluates avoiding stroke 
saturation through addition of a second low-order, high-stroke DM in tandem with the MEMS 
DM. 

2. Background 

MEMS devices were chosen as the high-order wavefront corrector for GPI because of their 
cost-effectiveness, high-actuator density, and agreeable performance. The mirrors can be batch- 
fabricated using a silicon foundry approach similar to that used to make integrated circuits, for 
a cost an order of magnitude cheaper per actuator than that of traditional piezoelectric DMs. 
MEMS technology is scalable to the actuator count required for ExAO high-contrast imaging. 

A series of prototype MEMS devices (two 144-element, ten 1024-element, and two 4096- 
element DMs) has been thoroughly tested at the University of California at Santa Cruz's Lab- 
oratory for Adaptive Optics (LAO). At the sub-nanometer level, hysteresis is negligible lH and 
MEMS devices are temporally stable and precisely positionablef9l. A MEMS DM has been 
flattened to 0.54 nm rms in the control band lflOl and has achieved 10^^ contrast at the far 
field ifTTI . An AO system using a 144-actuator MEMS DM has been successfully tested on- 
sky at Lick Observatorv 1121 \T3l . Figure [T] shows the LAO's first 4096-actuator DM, called an 
"engineering-grade" device because it is useable for testing and characterizing but is not fully 
functional across the mirror as required for GPI's "science-grade" MEMS. 

Residual low-frequency wavefront error in the GPI optical train will scatter star light into the 
region of the focal plane that should be kept dark in order to detect planets. This is often referred 
to as the "dark hole" region, and its size is set by the smallest-spatial-frequency control band, 
i.e. the spacing of the MEMS actuators in the Fourier domain lfT4l . Actuators on the deformable 
mirror that are broken, dead, or stuck scatter light into the dark hole by virtue of their influence 
functions producing a wide PSF in the far-field image. Saturated actuators are problematic for 
high-contrast imaging because, optically, insufficient stroke has the same effect as a broken, 
dead, or stuck actuator. 




Fig. 1. 4096-element MEMS deformable mirror in tlie Laboratory for Adaptive Optics. 
Center square is the reflective surface. Tlie DM is packaged on a ceramic chip carrier and 
enclosed within a refillable dry-nitrogen-gas chamber behind a pressure window. 

The following example (Fig. |2]i shows the effects of actuator saturation. A high-Strehl AO 
correction is effected by high-pass-filtering a Kolmogorov phase screen (left top, linear scale), 
with the corresponding image at the far field (left bottom, log scale). The far-field image is 
generated by the Fraunhofer approximation, propagating light from the phase-aberrated pupil 
plane (the DM) with an apodized transmission function ifTSlfTSl to the focal plane. If the phase 
measured at the MEMS plane is (j) and the apodizer function is A, then the intensity / at the 
far-field is simulated by calculating 



where A is the wavelength of the monochromatic light and FT means to take the Fourier trans- 
form. The second case (center) is an exaggerated illustration of actuator saturation. A low -pass- 
filtered Kolmogorov phase screen is clipped such that all values above 2 fim are set to 2 jxm. 
The residuals are displayed in linear scale (center top), and below that is the corresponding 
simulated far-field image in log scale. At right, radial averages of the two far-field images show 
the dramatic contrast achieved in the dark hole at low spatial frequencies when there is a high- 
Strehl AO correction (solid line). Therefore, stroke saturation (dashed line) cannot be tolerated 
since it obscures the dark hole and any planets located there. 

The statistics of Kolmogorov turbulence determine the stroke required of the wavefront cor- 
rector, whereas physical constraints on the DM determine its stroke capabilities. Variables that 
affect stroke capabilities include actuator design, size and spacing of actuators (i.e. pitch), and 
maximum voltage applied to the actuator. Due to the small (340-400-jLtm) pitch and the contin- 
uous facesheet on the MEMS, bending stresses limit high-spatial-frequency stroke on current 
devices to 1.2 /xm surface. The capacitor gap limits low-spatial-frequency stroke to 3.4 jj.m 
surface. Furthermore, stroke may also be drawn on to remove inherent curvature of the device 
and for tip/tilt residuals. 

Previous LAO articles reported on MEMS stroke as a function of spatial frequencvll9l [T7l[T8l . 
Three MEMS devices have been tested and the GPI MEMS has been specified for stroke at low 
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Fig. 2. Illustration of the consequences of stroke saturation. Left top: High-pass-filtered 
Kolmogorov phase screen, linear scale, and its simulated far-field image (left bottom), log 
scale. Center top: Exaggerated simulation of saturation with clipped Kolmogorov phase 
screen, linear scale, and its simulated far-field image (center bottom), log scale. Right: 
Radial averages of the simulated far-field images. Solid line: High-Strehl AO correction. 
Dashed line: Stroke saturated above 2 /im. 



and high spatial frequencies; Table[T]compiles the results. Norton et a/. lfTSl give more details of 
the "engineering-grade" 4k-MEMS characterization. The GPI "science-grade" 4k-MEMS has 
not yet been delivered, but its specifications are given. 
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Table 1 . Stroke measurements for a variety of MEMS DMs in the Laboratory for Adaptive 
Optics jT?). 'W107#X at a maximum voltage of 200 V is the same MEMS device used for 
the Kolmogorov saturation experiments described in Section|4] 



To measure MEMS stroke at a range of spatial frequencies, one-dimensional sinusoids of 
varying period were applied to two 1024-actuatorDMs, and the resulting stroke was measured 
in each easeful El- The results are plotted in Fig.[3]for the devices WI07#X (pitch 340 jxm 
and tested on the voltage range 0-160 V) and W95#39 (pitch 400 /im and tested on the voltage 
range 0-200 V). MEMS stroke decreases with increasing spatial frequency, whereas a hori- 
zontal line on Fig. [3] would be produced by both a segmented mirror and a DM with actuators 
having narrow sine influence functions (with the first zero occurring at a spacing of one actua- 
tor away). This departure from the narrow-sinc case indicates that the broad MEMS influence 
functions lfT9l have the effect of reducing the stroke as spatial frequency increases. 

Rather than plotting a single line with each data series on Fig. [3] two line segments joined 



by a "knee" fit each measurement series better The knee in both cases is located around 9 
cycles per aperture, a 3. 6-actuator period, while the influence function measured for a 340-/im- 
pitch MEMS falls to 4% at 2 actuators awav llT9l . about half this period. This implies that the 
high-spatial-frequency stroke is controlled locally by neighboring-actuators' influence func- 
tions, whereas the low-spatial-frequency stroke is in a distinct whole-MEMS regime. 
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Fig. 3. LAO testing of stroke as a function of spatial frequency for two MEMS deformable 
mirrors. Diamonds; Device W107#X, pitch 340 ^m, voltage range 0-160 V. Plus-signs: 
Device W95#39, pitch 400 flm, voltage range 0-200 V. The "knee" is located around 9 
cycles per aperture for both devices. Note that a DM having sine influence functions that 
fall to zero one actuator away would produce a horizontal line on this plot, as would a 
segmented mirror. 



3. Theory 

In a woofer-tweeter architecture, two DMs are arranged in series and both are conjugate to the 
same optical plane, usually that of the telescope pupil. The low-actuator-count, high-stroke DM 
is termed the "woofer," and the high-actuator-count, low-stroke DM is the "tweeter," a MEMS 
in the case of GPL As a consequence of the Kolmogorov turbulence spectrum, atmospheric tur- 
bulence strength decreases with spatial frequency 1201 . making the woofer-tweeter architecture 
feasible. The woofer DM corrects the higher-stroke, lower-spatial-frequency aberrations while 
the tweeter corrects the lower-stroke, higher-spatial-frequency aberrations. 

From Kolmogorov turbulence theory one can calculate the expected stroke requirement for a 
deformable mirror at a given telescope[2r|. At Gemini South, the planned location of GPI, the 
primary is D = 8 m and a typical Fried parameter is ro = 14 cm. The power spectral density <t> 
of Kolmogorov turbulence goes as <t>(K") 0= K"^^/^, where K is the wave number The formula 



for fitting error in radians Op as a function of ro and the actuator spacing d is 




where OCf is a constant that depends on the influence function for the deformable mirror in 
question ll22ll23l . 

For a continuous facesheet DM, Uf = 0.14 for circular segments with piston and tip/tilt 
actuation capabilities ll24l . Let us take this value for Uf, put the primary diameter D — d — %m, 
ro = 14 cm, and X = 500 nm. To calculate the total stroke requirement, we convert Equation|2] 
from square radians to microns surface, and multiply by 5 to obtain a ±5(7 correction: 

CTf = (±5) ^oF^^^ = ±2.17 ^m surface = 4.3 /im surface P-V. (3) 

The factor of 1 /2 converts from phase to surface and A /{2n) converts from radians to microns. 
However, 4.3 jim surface is the stroke requirement at the center of the pupil, and due to factors 
including non-stationarity of phase 1251 (which adds 20% at the edges of the pupil), inherent 
curvature of the MEMS (on the order of a few jim for current devices), and tip/tilt residuals, up 
to ^ 10 /xm stroke is desired 1261. From Table[T]we see that the 4k-MEMS devices do not get 
as much as 10 /xm stroke, so it was determined that a woofer DM was needed. 

Equation |2] can also be used to calculate the residuals after applying a woofer, which gives 
the stroke requirement for the tweeter in microns surface: 

where dn-^uf is the woofer actuator spacing. In this case the ap parameter becomes 0.28 — 
appropriate for the influence function of a continuous facesheet DM||24|. 

Varying d„„of in Equation |4] and multiplying the final result by ±3-5 to get the 3-5-sigma 
stroke requirement, the tweeter stroke required for each woofer degree of freedom is plotted in 
Fig.|4] Note that the stroke requirement quoted above of 10 /im is for a 8-m pupil and ro = 14 cm 
while Fig. |4]is for a 5-m pupil and ro = 15 cm. For the parameters in Fig. |4] using a woofer 
mitigates the tweeter stroke requirement to a more manageable ^1/im surface. Total stroke 
required rises as D6 with primary diameter (i.e. Equation|3]l, which increases the requirement 
on the woofer but not the tweeter. Rather, tweeter stroke is a function of the woofer pitch, 
not the primary diameter Therefore, experiments simulating a 5-m pupil are still valid for the 
8-m pupil case in regards to the tweeter. Laboratory experiments follow to verify empirically 
the ability of the MEMS to form Kolmogorov atmosphere shapes, using a woofer to avoid 
saturation. 



4. Experiment 

We tested MEMS device W107#X for stroke saturation when being deformed to correct for 
Kolmogorov turbulence. W107#X has a 340-jLim pitch and was applied with a maximum volt- 
age of 200 V. From Table [H we see its low-order stroke is 1.5 /xm surface and its high-order 
stroke is 0.2 jim surface. While this DM does not have as much stroke as the final GPl 4k- 
MEMS will, the experiments are useful in determining whether inter-actuator or peak-to-valley 
stroke causes more saturation events. 

Device W107#X has 1024 actuators (32 across) rather than the 44-actuators-across GPI will 
use. Thus, because the high spatial frequencies on the MEMS were of more interest than the 




Fig. 4. Tweeter stroke required as a function of woofer degrees of freedom for ro = 15 cm, 
based on Equation |4] 



low, the MEMS pitch with respect to the primary was held fixed at the GPI value (18 cm), mean- 
ing the size of the pupil in this experiment corresponds to a 5-m telescope. The comparative 
mapping between the MEMS actuators and the telescope pupils is shown in Fig. |5] superim- 
posed over a Kolmogorov phase screen. The entire square image contains 44x44 actuators, 
the outer inscribed circle delineates an 8-m pupil (GPI), the inner black square marks 32x32 
actuators, and the inner inscribed circle delineates a 5-m pupil (this laboratory experiment). 

Measurements were made using the phase-shifting diffraction interferometert27il , a sub- 
nanometer-absolute-accuracy measuring instrument situated on the ExAO testbed ll28l l29ll at 
the LAO. A spherical wave (f/220) converges from a beam size of 12 mm at the MEMS plane 
to the surface of a pinhole aligner, where it is combined with a reference beam. The reference 
beam emerges from a fiber in spherical waves, emitted through the l.S-fim pinhole to interfere 
with the reflected test beam. This reference wave has an excellent wavefront quality (A /500) 
with which to compare the test beam. The resulting diverging spherical wave fringes are imaged 
on a CCD camera. The interferogram is numerically back-propagated using a Huygens method 
to obtain phase and amplitude data at the MEMS plane. 

The free parameters in the saturation experiment were the Fried coherence size, ro, and the 
woofer pitch, d^voof- The procedure was as follows: An atmospheric turbulence screen with the 
specified ro was generated in software by enforcing Kolmogorov turbulence statistics upon an 
array of random numbers. Compensation of the atmosphere screen by an ideal woofer DM with 
the specified d„>„o/ was simulated by removing all the power at low spatial frequencies below 
the appropriate cutoff frequency in the Fourier domain. The resultant phase screen contained 
only the higher spatial frequencies. The MEMS was then "closed to" this turbulence screen 
in closed-loop by using the screen as the reference. That is, rather than taking a flat shape as 
the target of the closed-loop algorithm, the target for convergence was the phase screen. After 




Fig. 5. To-scale mapping of 8m Gemini pupil (outer gray circle), 5m testbed pupil (inner 
gray circle), 44x44-actuator array for GPI (full grid), and 32x32-actuator array of MEMS 
in testbed (black). Fine white grid delineates individual MEMS actuators. 



running 20 closed-loop iterations, the iteration with the lowest rms wavefront error between the 
target phase screen and the measurement was taken for analysis, along with the corresponding 
voltage array on the MEMS. Convergence did not generally improve after 15 iterations. 

Fig.|6]is a histogram of the stroke inputs at individual actuators' locations. Three curves are 
plotted, showing the input Kolmogorov phase screens for ro = 15 cm cases with no woofer 
(piston and tip/tilt removal only), 100-cm-pitch woofer, and 50-cm-pitch woofer. There were 
ten screens of each case for these parameters. This plot shows the stroke range applied to the 
individual actuators to make the input Kolmogorov phase screens. The woofer dramatically 
reduces the phase variance input to the tweeter. 

5. Results 

Taking the best closed-loop iteration in each case, any phase- wrapping (an interferometry ar- 
tifact) was removed by masking out the bad data points. Additionally masked out was the 3x3 
region centered on a dead actuator located at coordinate [22, 8] (the origin is in the lower-left- 
hand corner and counting starts with 0). This was to avoid its being counted as "saturated" in 
each frame. 

The corresponding array of voltages applied to the MEMS was used to determine saturation 
by the criteria that any actuator commanded to the minimum (0 V) or maximum (200 V) was 
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Fig. 6. Histogram of individual actuator phases input for the ro = 15 cm experiment. 
Dashed line: No woofer (piston and tip/tilt removal only); Dashed-dotted line: 100-cm- 
pitch woofer correction; Solid line: 50-cm-pitch woofer correction. 

declared "saturated." Table |2] lists the set of parameters explored and gives the number of ac- 
tuators saturated. The number of actuators tested varies due to different numbers of trials and 
different numbers of actuators remaining within the unmasked region after discarding phase- 
wrapped data. For ro = 10 cm and d„>„o/ = 31 cm, zero actuators were saturated out of the 4852 
tested. Therefore, the probability of saturation in that case was calculated using the binomial 
theorem and taking the outside-3-(7 probability for the upper limit. Figure [T] plots the satura- 
tion frequency as a function of woofer pitch for ro = 10-15 cm. Error bars are standard error 
{o /\/N). a power law is fit to each series. 

Actuator voltage maps give a qualitative picture of close saturation encounters for the MEMS 
in closing to a Kolmogorov phase screen. Figure [H] shows three typical actuator voltage maps. 
Saturated actuators are identified with a white cross (saturated low) or a black diamond (sat- 
urated high). The first case is a typical iteration for iq = 10 cm and no woofer (piston and 
tip/tilt removal only). The mean number of saturated actuators for these parameters was 16 per 
trial, and this particular iteration has three actuators saturated high and eleven saturated low. 
The second case is a typical iteration for ro = 10 cm and d„,„o/ =134 cm. The mean number 
of saturated actuators for these parameters was two per trial, and this particular iteration has 
two actuators saturated low and none saturated high. The third case is a typical iteration for 
ro = 10 cm and d„(,of = 62 cm. The mean number of saturated actuators for these parameters 
was one per trial, and this particular iteration has one actuator saturated low and none saturated 
high. 

Besides using these actuator voltage maps to display the number of saturated actuators, we 
also notice something qualitative yet important with these plots: saturated actuators occur either 
singly or else in clumps. They do not, however, occur in pairs of one up and one down: in the 
clumpy patches of saturated actuators, all members of the saturated clump are either high or low. 
In all, 140 trials total were done (57 046 actuators were tested after masking out bad data), and 
there are no instances of two adjacent saturated actuators in which one was saturated high and 
the other low. Thus, the most extreme case of inter-actuator saturation at opposing positions 
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Fig. 7. Saturation frequency as a function of woofer pitcfi. Data for experiments at ro = 10- 
15 cm are represented by asterisks and diamonds, respectively, while a power-law fit to each 
experiment is represented by a dotted and a dashed line, respectively. Solid and dashed- 
dotted curves are predictions for ro = 10-15 cm, respectively, calculated with Equation|4] 
(see details in Section[6}. The fitting-error constant af used to derive the predictions was 
varied smoothly from 0.14 (piston/tip/tilt only) to 0.28 (continuous facesheet DM) over the 
range d„.„oy/D > 0.25. This is due to mapping arguments: at the largest d^^.ogffD ratios, the 
pitch of one woofer actuator is a good fraction of the primary diameter, and thus whether 
the woofer actuators overlap or align with the pupil edge becomes important in determining 
the fitting error. This choice (of 4-woofer-actuators-across as the limit for where the DM 
becomes like a continuous facesheet) affects the steepness of the predicted curves but not 
their endpoints, and we note that only one of our experimental data points (290-cm-pitch 
woofer for ro=15 cm) is well inside this regime where Uf is approximated. We see that the 
saturation incidence if narrow sine influence functions are assumed should have fallen off 
more steeply than measured. Thus decreasing the woofer pitch did not reduce saturation 
occurrence as much as expected. Moderate to high spatial frequencies are saturating more 
on the MEMS than expected, and a simple analytical calculation underpredicts the stroke 
requirement. 
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Table 2. Stroke saturation results: Fraction of actuators saturated for each set of (rg, d^„jgf ) 
parameters, after removing piston and tip/tilt. The number of actuators tested in each case is 
equal to the total number of MEMS actuators within the pupil, summed over all trials. This 
quantity varies according to two effects: there is a non-uniform number of trials conducted 
for each set of parameters, and there is a non-uniform number of actuators analyzed inside 
the pupil for each trial. The pupil nominally encircles ~ 600 actuators, but actuators are 
masked out and not counted if they are in regions where phase-wrapping occurs in the 
interferometric measurement. *At the largest woofer-pitch end, no woofer was used — only 
piston and tip/tilt were removed over a 5-m aperture. We therefore adopted "500 cm" as 
the abscissa for these data to be able to plot the corresponding ordinate from column 4 onto 
Fig.|7] 'Upper limit. 

is not occurring. However, we will see below that moderate to high spatial frequencies were 
problematic in another sense. 

Figure |9] shows the peak-to-valley stroke histogram for various woofer pitches, as noted in 
the legend. The trials without a woofer (piston and tip/tilt removal only) had a peak-to-valley 
stroke ranging from 0.8 to 1.7 fim surface. Using a woofer reduces the peak-to-valley stroke 
required of the MEMS by a factor of 3-4, down to 0.2 to 0.5 /xm surface. Note that a maximum 
of 1.5 jim surface stroke was measured for this MEMS, device W107#X (Table[TJ, in the low- 
frequency stroke test. Obtaining more stroke, up to 1.7 jim surface, with a Kolmogorov phase 
screen is attributed to the curvature of the inherent shape of the MEMS, which is 450 nm surface 
for this device. 

Figure[TO]shows the inter-actuator stroke histogram for the same cases. Inter-actuator stroke 
is calculated by differencing the measured phase between each pair of neighboring actuators. 
Neighbors to a particular actuator are horizontally, vertically, or diagonally adjacent to that 
actuator The woofer reduces the required inter-actuator stroke by a factor of 2-3, a smaller 
improvement than the peak-to-valley reduction. This is as expected because the woofer actuator 
spacing remains much larger than the inter-actuator MEMS spacing (340 jim), so the inter- 
actuator stroke is at a higher spatial frequency than the woofer can control. In addition, a high- 
amplitude low-frequency aberration will create a slope between two neighboring actuators, 
introducing further power at the inter-actuator level. 

6. Analysis 

In order to be able to scale these results to other configurations than our ro and d,„,o/ ranges on 
a 5m primary, we design a model and test it by using it to predict the measured saturation oc- 
currence. A histogram of Kolmogorov-turbulence phase (i.e. Fig.|6]l follows Gaussian statistics, 
which can be plotted from theory by calculating the standard deviation with Equation |4] For 
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Fig. 8. Maps of MEMS actuator voltages, with saturated actuators marked as follows: ac- 
tuators saturated low are marked with white crosses, whereas actuators saturated high are 
marked with black diamonds. Typical results are shown. Left: vq = 10 cm, no woofer (pis- 
ton and tip/tilt removal only). Three actuators were saturated high and 11 were saturated 
low. Center: ro = 10 cm, d„.„oy =134 cm. Two actuators were saturated low (none were 
saturated high). Right: ro = 10 cm, A„oof = 62 cm. One actuator was saturated low (none 
were saturated high). 

a DM having sine influence funetions with the first zero located one actuator away, the stroke 
of the DM is the same at all spatial frequencies. Suppose the stroke of such a sinc-influence- 
function DM is 1.5 /im surface. Plot a histogram of Kolmogorov turbulence, and the phase 
within +1.5 /im and —1.5 /xm are within the capabilities of this sinc-influence-function DM to 
correct, whereas the area outside the ±1.5 fim cut-off is proportional to the fraction of actuators 
that saturate. 

Table [T] and Fig. [3] show, however, that the MEMS influence functions are not ideal. There- 
fore, instead of using a constant 1.5 jxm as the stroke cut-off, we can take the stroke at each 
spatial frequency, scaling it by a constant multiplier to go from the measurements at 160V to 
200V. Table[3]lists the predicted number of actuators saturated for each set of parameters using 
these methods of finding the area under the Gaussian distribution beyond the cut-off stroke. 
Two cut-off strokes are used for comparison, the ideal DM case and the MEMS W107#X case 
with stroke varying as a function of spatial frequency. The measurements from Table |2] are 
repeated for comparison. 

Saturation was as frequent as expected for cases with larger woofer pitches, but was more 
frequent than expected for cases with smaller woofer pitches. The moderate to high spatial 
frequencies were predicted to cause less saturation on the 1 .5-/im-stroke MEMS than occurred. 
Therefore, since the model did not fit the data in this case, the theory could not be used to 
extrapolate to other configurations. However, the GPI 4k-MEMS is specified to have 5 times 
more high-order stroke than the MEMS used in this experiment, which will further reduce 
stroke saturation. 

Surprisingly, even the second attempt at prediction (using a variable stroke cut-off for each 
particular spatial frequency) did not match the experiments at the moderate to high spatial 
frequencies. This can be explained by looking again at Fig.[3]and seeing how stroke continues 
to fall off, albeit less steeply past the knee, no matter what spatial frequency is taken. If the 
stroke beyond a certain point were uniform, we would be in the ideal-influence-function regime 
for those spatial frequencies. The prediction, by using a single stroke value (at each spatial 
frequency) for the cut-off of the Gaussian distribution, is still assuming uncorrelated narrow- 
sinc-function actuators, whereas on the MEMS there are always higher spatial frequencies that 
limit the stroke. 
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Fig. 9. Histogram of peak-to-valley (i.e. low-spatial-frequency) stroke measured for all 
MEMS actuators in the ry = 15 cm experiment. Dashed line: No woofer (piston and tip/tilt 
removal only); Dashed-dotted line: 100-cm-pitch woofer correction; Solid line: 50-cm- 
pitch woofer correction. The woofer reduces the peak-to- valley stroke on the MEMS by a 
factor of 3-4. 




Fig. 10. Histogram of inter-actuator (i.e. high-spatial-frequency) stroke measured for all 
MEMS actuators in the ro = 15 cm experiment. Dashed line: No woofer (piston and tip/tilt 
removal only); Dashed-dotted line: 100-cm-pitch woofer correction; Solid line: 50-cm- 
pitch woofer correction. The woofer reduces the inter-actuator stroke on the MEMS by a 
factor of 2-3. 
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Table 3. Table of measured versus predicted stroke saturation occurrence. Predictions were 
made by using Equation |4] as the standard deviation to plot a Gaussian distribution of all 
the probable Kolmogorov phase values, taking the area under the curve outside the cut-off 
stroke to be proportional to the fraction of actuators saturated. The failure of the predicted 
values to match the measured values, except at the lowest spatial frequencies, indicates 
that the decreasing stroke with increasing spatial frequency of W107#X is producing more 
saturation than expected for a segmented DM or a DM with sine influence functions. *At 
the largest woofer-pitch end, no woofer was used — only piston and tip/tilt were removed 
over a 5-m aperture. We therefore adopted "500 cm" as the abscissa for these data to be 
able to plot the corresponding ordinate from columns 3-5 onto Fig.|7] ' Upper limit. 

7. Discussion 

In systems where saturation cannot be tolerated, the wavefront correctors must have sufficient 
stroke to correct the amplitude of the phase aberrations. If deformable mirror stroke is lim- 
ited, stroke can be conserved for dynamical corrections by instead using static correctors to 
remove inherent curvature to the DM itself and for correcting static aberrations in the opti- 
cal path. Similarly, the unpowered shape of the DM itself should be as flat as possible so as 
to conserve stroke for correcting turbulence. The stroke requirement increases with power in 
the Kolmogorov spectrum, so better astronomical seeing conditions will give less saturation. 
Accordingly, the astronomical site and observing conditions could be favorably selected to 
minimize stroke saturation. 

In designing a deformable mirror, factors that reduce stroke saturation include decreasing the 
thickness of the facesheet, increasing the spacing (pitch) between actuators, increasing the op- 
erating voltage, and modifying the mechanical design of the actuator (varying parameters such 
as the capacitor gap, spring constant, flexural rigidity of the plate, or using leverage bending). 
These are varied primarily at the research and development phase. 

In this work we assume that piston and tip/tilt modes are corrected by other devices, and that 
the stroke of the woofer DM is high enough to correct the frequencies below the cutoff; there- 
fore, we are only concerned with avoiding saturation on the tweeter DM. In a dual-DM system, 
factors that reduce tweeter saturation via sending saturable stroke to a second low-order woofer 
DM include: employing a woofer with more actuators across the pupil (the subject explored in 
this paper); and carefully controlling the frequency distribution sent to each wavefront corrector 
in the control loop 1261 l30l . 

Tweeter saturation occurrence can be predicted by analyzing simulations of correcting turbu- 
lence screens with model MEMS mirrors llSTl [8l [32l [33l . Our results show the analytic fitting- 



error model based on the turbulence power spectrum underpredicts incidence of stroke satura- 
tion. Our original goal in this work was to measure whether stroke saturation would be prob- 
lematic for GPL In order to assure the higher-stroke MEMS mirrors (now being fabricated) are 
adequate for GPI, we plan to measure saturation statistics with similar testing on "engineering- 
grade" versions of the 4k-MEMS. 

8. Conclusions 

In this experiment we applied Kolmogorov-turbulence atmospheric phase screens to a MEMS 
deformable mirror The ro of the phase screens ranged from 10-15 cm, and the pitch of the 
woofer used to remove the low spatial frequencies ranged from 8-16 actuators across an 8-m 
primary. 

The MEMS when solitary suffered saturation ^4% of the time. Using a woofer DM reduces 
dramatically the amount of phase sent to the MEMS, and thus mitigated MEMS saturation 
occurrence to a fraction of a percent. 

The woofer did not mitigate mid-to-high-spatial-frequency stroke as much as expected. No 
neighboring actuators were saturated at opposing positions, meaning the most extreme case of 
inter-actuator saturation did not occur 

However, moderate to high spatial frequencies did saturate more often than predicted based 
on a Gaussian distribution of phase and somewhat-idealized influence functions. This implies 
that correlations in actuators, i.e. high-spatial-frequency stroke limits, are significant. Attempts 
to derive the expected saturation using a simple analytical model underpredicted the stroke 
requirements, showing that empirical studies are important. 
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